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Abstract: We study various aspects of the matrix models calculating free energies and Wilson 
loop observables in supersymmetric Chern-Simons-matter theories on the three-sphere. We first 
develop techniques to extract strong coupling results directly from the spectral curve describing 
the large N master field. We show that the strong coupling limit of the gauge theory corresponds 
to the so-called tropical limit of the spectral curve. In this limit, the curve degenerates to a planar 
graph, and matrix model calculations reduce to elementary line integrals along the graph. As 
an important physical application of these tropical techniques, we study M = 3 theories with 
fundamental matter, both in the quenched and in the unquenched regimes. We calculate the 
exact spectral curve in the Veneziano limit, and we evaluate the planar free energy and Wilson 
loop observables at strong coupling by using tropical geometry. The results are in agreement 
with the predictions of the AdS duals involving tri-Sasakian manifolds. 
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1. Introduction 

Recently, the localization methods of [1] were extended to supersymmetric Chern-Simons-matter 
theories in [2, 3, 4]. These methods provide matrix model representations for the Euclidean 
partition function of these theories on the three-sphere, as well as of BPS Wilson loop VEVs. In 
the case of the ABJM theory constructed in [5, 6], the resulting matrix models were solved at all 
orders in the 1/N expansion [7, 8] by using previous results [9, 10, 11] on Chern-Simons matrix 
models [12] and their large N topological string duals. This exact solution made possible the 
calculation of exact interpolating functions in ABJM theory, testing in this way various strong 
coupling predictions based on its AdS4 dual. In particular, the A 3 / 2 behavior of the number of 
degrees of freedom of the M2 brane theory [13] was derived in [8] by calculating the exact planar 
free energy on the three-sphere. 

The Chern-Simons-matter quiver theory of ABJM has been generalized in many directions, 
by considering for example more complicated quivers or by adding flavor, i.e. matter multiplets in 
the fundamental representation. The resulting theories have conjectural AdS duals which lead to 
predictions for the strong coupling behavior of various observables. If the Chern-Simons-matter 
theories have enough supersymmetry, the methods of [2] provide matrix model representations 
for the partition function on the three-sphere and for Wilson loop VEVs. In general, these matrix 
models are quite complicated, and solving them for all values of the coupling, even in the planar 
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limit, is a difficult task. In this respect, the ABJM theory is rather exceptional since it can 
be related to a well-known topological string theory, and this is what made possible to obtain 
very explicit interpolating functions [7, 8]. In some more general cases, like the one considered 
in [14], one can write down the planar resolvent in terms of contour integrals, but the resulting 
expressions are not explicit enough to extract in a straightforward way the physical information 
on the observables. 

As a matter of fact, finding an exact solution to the matrix model at all values of the coupling 
might be a Pyrrhic victory: if we just want to test the AdS predictions for a given Chern-Simons- 
matter theory, we simply need the leading result at strong coupling. Therefore, we would like to 
have a method to perform calculations in the matrix model directly in the strong coupling limit, 
in the same way that we can use matrix model perturbation theory to obtain results in the weak 
coupling limit. 

The first goal of this paper is to develop techniques to perform such strong coupling compu- 
tations. The basic idea is that, as emphasized in [8], the strong coupling behavior of the ABJM 
theory corresponds to the limit of large complex structure for the spectral curve describing the 
large N distribution of eigenvalues. It turns out that this limit has a number of simplifying 
features that have been much studied recently in the mathematical literature: it is the so-called 
tropical limit of the curve. In this limit, the spectral curve degenerates to a two-dimensional 
graph, and the calculation of periods reduces to the computation of line integrals along this 
graph. Many recent results in tropical geometry can in fact be adapted to our purposes, and we 
obtain in this way simple methods to evaluate the strong coupling limit of observables in ABJM 
theory. 

The main reason to develop techniques based on tropical geometry is of course to be able 
to go beyond ABJM theory and obtain strong coupling results for more complicated Chern- 
Simons-matter theories from their matrix models. In this paper we will focus on the theories 
with fundamental matter introduced in [15, 16, 17]. These theories have only M = 3 super- 
symmetry, and they are simple toy models for the incorporation of flavors in the AdS/CFT 
correspondence. As in other theories with fundamental matter, one can study them in the 
quenched approximation, or one can consider the full unquenched regime, in which matter loops 
are not suppressed. We show that both approaches can be incorporated in the matrix models 
of [2] . The quenched approximation leads to a systematic expansion in powers of the number of 
flavors, where each term can be computed as a correlator in ABJM theory. In order to study 
the unquenched theory we find the exact planar resolvent of the N = 3 theory in the Veneziano 
limit [18]. In both cases, quenched and unquenched, we can use tropical techniques to calculate 
observables at strong coupling, testing in this way the AdS predictions for these theories made 
in [15, 16, 17]. In particular, we reconstruct from the matrix model calculation the full volume 
of the tri-Sasaki Einstein manifold appearing in their M-theory dual. 

The organization of this paper is as follows. In section 2 we review the relevant Chern- 
Simons-matter theories, their matrix model description, as well as their AdS duals. In section 
3 we review the exact solution of the ABJM theory matrix model worked out in [7, 8], and we 
show how to perform direct strong coupling calculations by using the tropical limit of the exact 
solution. In section 4 we start our study of the N = 3 Chern-Simons-matter theories with flavors 
and we present the matrix model analogue of the quenched expansion. We also evaluate the first 
terms in this expansion at strong coupling by developing techniques based on tropical geometry. 
In section 5 we solve the N = 3 theory in the Veneziano limit, we test the exact planar resolvent 
at weak coupling, and we calculate the free energy and Wilson loop VEVs at strong coupling, 
recovering in this way the AdS predictions. Finally, in section 6 we state our conclusions and 
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open problems. 

When the first version of this paper was being typed, the very interesting paper [19] appeared, 
which has some overlap with our results. One of the goals of [19] is also to obtain results directly 
at strong coupling in the matrix model, and they develop a simple and powerful technique to do 
that based on a careful analysis of the eigenvalue density. Their main conclusion, that the density 
of eigenvalues is piecewise linear in this regime, corresponds precisely to the tropical limit of the 
spectral curves considered in this paper. In the Appendix we use the methods of [19] to re-derive 
the free energy of the N = 3 theory with fundamental matter, and we clarify the relationship 
between their approach and the tropical techniques of this paper. 



2. Supersymmetric Chern— Simons— matter theories and matrix models 
2.1 ABJM theory and matrix models 

The ABJM theory [5, 6] is a supersymmetric quiver Chern-Simons-matter theory in three di- 
mensions with gauge group U(N\)k x U{N-2)-k and M = 6 supersymmetry. The Chern-Simons 
actions have couplings k and —k, respectively. The theory contains in addition four hypermulti- 
plets: there are two bifundamental chiral fields (A±, A 2 ) in the representation (N\, N2), and two 
bifundamental chiral fields {B\,B2) in the representation (N±, N 2 ). 

The ABJM matrix model, obtained in [2] , gives an explicit integral expression for the parti- 
tion function of the ABJM theory on S 3 , as well as for Wilson loop VEVs. This matrix model is 
defined by the partition function 

Zabjm(Ai, N 2 ,g s ) 

(2.1) 

where the coupling g s is related to the Chern-Simons coupling k of the ABJM theory as 

9s = "f. (2.2) 

The normalization of this matrix integral has been explained in [8] . It was shown in [7] that the 
ABJM matrix model is closely related to the L(2, 1) lens space matrix model introduced in [12] 
and further studied in [9, 10, 20]. This matrix model is defined by the partition function 



w*. = wWn t n t n H (^)) 2 H 



n(2cosh(^p)) e-^G^+EW). 

(2.3) 



x 

';3 



The relation between the partition functions is simply 

Z A BJM(Aq,Ar 2 ,<7 s ) = Z £(2)1) (Wi,-JV2,<fc). (2.4) 

Since the large N expansion of the free energy gives a sequence of analytic functions of N±, N2, 
once these functions are known in one model, they can be obtained in the other by the trivial 
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change of sign N2 — > — N2. Equivalently, it can be shown [7, 4] that the ABJM matrix model is 
the supergroup version of the lens space matrix model of [12, 9]. 

Natural observables in this theory are supersymmetric Wilson loop operators. The 1/6 BPS 
Wilson loop has been constructed in [21, 22, 23] and it is labelled by a representation R of 
U(N\). One of the main results of [2] is that the VEV of this Wilson loop in ABJM theory can 
be obtained by calculating the normalized VEV of the matrix e w in the matrix model (2.1), i.e., 

(W 1 R /6 )=g s (Tv R (e^)) ABm , (2.5) 

1/2 

A 1/2 BPS loop W% was constructed in [4] , where 1Z is a representation of the supergroup 
U(Ni\N2). In [4] it was also shown that it localizes to the matrix model correlator in the ABJM 
matrix model 

(W% 2 ) = g s (Sti n U) ABm , (2.6) 



U=[ e "L)- (2-7) 



where U is the diagonal matrix 

-e^ 

2.2 Flavored theory and its gravity dual 

It is possible to flavor the ABJM theory by adding matter hypermultiplets in the fundamental 
representation [15, 16, 17]. More precisely, one adds multiplets (Qi,Qi), with i = 1,2. The 

fields Qi, i = 1,2 are in the representations (A^,l) and (1,^), respectively, while Qi are in 
the conjugate representations (Ni, 1) and (1, A/" 2 ), respectively. This matter content breaks the 
M = 6 supersymmetry of the ABJM theory down to M = 3. Notice that the ABJM theory can 
be obtained, formally, as the limit 

Nf } -»• (2.8) 

of the flavored theory. We will denote by 

N f = N { f 1} + Nf ] (2.9) 

the total number of flavours. 

The partition function of this N = 3 theory can be also computed in terms of matrix models, 
as explained in [2]. The inclusion of extra matter hypermultiplets just leads to the insertion of 
determinant-type operators in the matrix integral (2.1): 



M g du 3 (2-h (^)) 2 (2sinh 



Nim j 11 27r n 2. n _. ( 2cogh (2.10) 

n (2 -sh ^y Nf) n (2 ^h |) ~ Nf) e-u^-^). 



x 

=1 



The large iV dual of the ABJM Chern-Simons-matter theory is given by type IIA string 
theory on AdS4 x CP 3 , which lifts to M-theory on AdS4 x § /Zfc [5, 6]. It was conjectured in 
[15, 16, 17] that, when iVi = A2 = N, the N = 3 theory with flavor has a type IIA large N dual 
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where Nf D6 branes wrap the MP 3 cycle inside CP 3 . This is the four- dimensional counterpart of 
the original construction of [24], which adds flavor to AdSs by wrapping D7 branes around an § 3 
inside S 5 . The flavored Af = 3 theory also describes N M2 branes probing an eight-dimensional 
hyperKahler cone Ais with Sp(2) holonomy. The base of this cone is a tri-Sasakian space X 7 . 
The space A4s is a particular member of a family of hyperKahler cones A4s(t) labeled by three 
natural numbers t = (ti,t2,ts). These cones can be constructed as hyperKahler quotients 

M 3 ///U(l), (2.11) 

where the U(l) action is characterized by the three charges t. The bases X 7 (t) of these cones give 
an infinite family of tri-Sasakian manifolds known as Eschenburg spaces, see [25] for a detailed 
study and references to the relevant literature. The dual to the N = 3 Chern-Simons-matter 
theory with a total number of Nf fundamentals has charges 

t = (Nf,N f ,k). (2.12) 

In the following, the eight-dimensional cone corresponding to this charge will be simply denoted 
by M 8 - 

At large N the above theory of N M2 branes is described by M-theory on the manifold 

AdS 4 x X 7 , (2.13) 

where X-j is the tri-Sasakian seven manifold corresponding to (2.12). This background is the 
eleven-dimensional lift of the type IIA background with Nf D6 branes. The eleven-dimensional 
metric and four-form flux are given by 



d«n = R\ 7 ( Td4ds 4 + ds 'xr I > 



G = -Rx 7 UAds 



(2.14) 



where WAdS 4 is the volume form with unit radius. The radius Rx 7 is determined by the flux 
quantization condition 

(2ir£ p fN= [ * u G = 6R X7 vo\(X 7 ). (2.15) 
Jc 7 

In this equation, £ p is the eleven-dimensional Planck length, C-j is a cycle enclosing the brane 
and homologous to X7, and the volume of X7 (with unit radius) is given by [25] 

vo.(X 7 ) = , ( , W , 

where 

V 1 + A*/2 k 

Newton's constant in four dimensions can be obtained by standard compactification of the 
Einstein-Hilbert action in eleven dimensions, 

1 k 1 / 2 iV 3 / 2 

-M. (2-18) 



16vrG 4 12V2tt 
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This conjectural large N dual of the N = 3 theory makes two quantitative predictions for the 
strongly coupled N = 3 Chern-Simons-matter theory with flavors. First of all, the free energy 
of the M = 3 theory on § 3 is given by the same expression that the one for the ABJM theory 
[26, 8] but with the new Newton's constant (2.18), 



7T 

2Ga 



N 3 / 2 k 1 ^ (fi) 



(2.19) 



Another quantity that we are interested in is the VEV of supersymmetric Wilson loops. As usual 
in the AdS/CFT correspondence, this can be calculated by evaluating the regularized area of a 
fundamental string in the type IIA reduction of the above M-theory background. The resulting 
geometry, which includes the full backreaction of the D6 branes, is a warped compactification 
and we have not performed such a calculation. However, it was pointed out in [16] that 



and we then expect 



- tl str 



1 



2-kN 1 / 2 



4N f + k ^jN f + 2k 



exp 



^jwy+zk 



(2.20) 



(2.21) 



for both the 1/2 and 1/6 BPS Wilson loops. This VEV incorporates the screening effect on 
Wilson loops due to unquenched flavor. Indeed, we see from (2.21) that, as the number of flavors 
grows, the exponent decreases. This might be interpreted as a conformal avatar of the screening 
effect. Also notice that, when Nf —> 0, one recovers in (2.21) the right value for the ABJM limit. 
Notice that, in the computation leading to (2.21), a possible contribution to the vev of strings 
ending on the D6 branes has not been taken into account. We will see however in this paper that 
a gauge theory computation at strong coupling leads to a result in agreement with (2.21), thus 
indicating that such contributions are absent or subleading. 

3. Strong coupling limit and tropical geometry 

The ABJM matrix model can be solved in the planar limit simply by borrowing the exact solution 
of the lens space matrix model worked out in [9, 10]. Let us now review some ingredients of this 
solution following the notations of [8]. 



7ri — B 
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Figure 1: Cuts in the z-plane and in the Z-plane. 
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At large N, the two sets of eigenvalues, fii, Vj in (2.3) condense around two cuts. The cut of 
the fj,i eigenvalues is centered around z = 0, while that of the vi eigenvalues is centered around 
z = 7ri. We will write these cuts as 

d = (-A,A), C 2 = (TTi-B,Tri + B), (3.1) 

in terms of the endpoints A, B. It is also useful to use the exponentiated variable 

Z = e z , (3.2) 

In the Z plane the cuts (3.1) get mapped to 

(l/o,a), (-l/b,-b), a = e A , b = e B , (3.3) 

which are centered around Z = 1, Z = — 1, respectively, see Fig. 1. We will use the same notation 
Ci,2 for the cuts in the Z plane. The total resolvent of the matrix model, co(z) is defined as [10] 

= 9s (Tr (§±|) ) = g. coth (^i) ) + ^ (£ tanh (^) ) ( 3 - 4 ) 

where U is given in (2.7). 

We will denote by ujq(z) the planar limit of the resolvent, which was found in explicit form 
in [10]. It reads, 

MZ)^- = log(^ [f(Z) - y/p{Z)-4(PZ*] ) ^, (3.5) 

where 

f(Z) = Z 2 - (Z + 1. (3.6) 
Here, /3 and C are "bare" parameters characterizing the resolvent. If we write 

a(Z) = f(Z) - 4/3 2 Z 2 = (Z -a)(Z- 1/a) [Z + b)(Z + 1/6) (3.7) 

we see that the parameters £, ft are related to the positions of the endpoints of the cuts as follows 

C = Ua + --b-l) , ^=]( a + 1 +b+l). (3.8) 



2 V a bj 4 V a 6, 

All the relevant planar quantities can be expressed in terms of period integrals of the one- form 
ujq{z)Az. The 't Hooft parameters 

U = 9sNi, * = 1,2, (3.9) 

are given by the period integrals 

U = — I to {z)dz, i = l,2. (3.10) 

The total 't Hooft parameter 

t = h + t 2 (3.11) 

is related to the modulus (3 as 

P = e*. (3.12) 
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The planar free energy Fq satisfies the equation 

dF dF If 
l= — - — - m t = --^ u (z)dz, (3.13) 

where the T> cycle encloses, in the Z plane, the interval between —1/6 and 1/a, see Fig. 1. 
Another quantity that can be computed as a period integral is the 1/6 BPS Wilson loop, since 

[7] 

The 1/2 BPS Wilson loop is then given by 

We finally recall that the solution to the original ABJM matrix model can be obtained by simply 
setting 

ti = 2vriAi, t 2 = -2vriA 2 , (3.16) 

where 

A, = f (3.17) 

are the 't Hooft parameters of the ABJM model. 

In order to make contact with the algebro-geometric formalism based on the spectral curve, 
notice that if we denote 

log Y = u {Z) (3.18) 
we find the following equation for an algebraic curve in C* x C*: 

Z 2 fofrrt 1 



y + _-^-MZ' Z + 1)=0, (3.19) 



where the parameters z±, z 2 are related to Q,j3 as 



C=-^, P = ^- (3-20) 



zi V z i 

The curve (3.19) is the mirror curve to the toric Calabi-Yau manifold known as local P 1 x P 1 , and 
the parameters z\ , z 2 become complex coordinates for the moduli space of the mirror Calabi-Yau. 

In order to make contact with the AdS dual, one has to calculate the gauge theory/matrix 
model quantities at strong coupling. In [7, 8] this was done essentially by computing exact 
interpolating functions at all couplings and then going to the strong coupling regime. However, 
the calculation of interpolating functions might become hard, specially in more complicated 
generalizations of ABJM theory like the theories with matter considered in [15, 16, 17]. In 
particular, one would like to have a computational framework to do calculations directly at 
strong coupling, without going through the determination of exact interpolating functions. We 
will now propose such a framework, and we will illustrate it by considering the ABJM theory. 

As explained in [8] , the strong coupling limit of the ABJM theory corresponds to the large 
radius limit of this Calabi-Yau moduli space. This is the limit where 

z u z 2 -> 0. (3.21) 
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It is well-known that in this regime there are two natural period integrals T\ , T 2 given by 



-Ti =logzi + J 1) {z 1 ,z 2 ), 
-T 2 = logz 2 +J 1 \z 1 ,z 2 ), 

where 



(3.22) 



J 1 Hz 1 ,z 2 ) = 2 £ r(i + + /) 2 ^ = 2zi + 222 + 3z ' + 12zi22 + 3z 2 + • • • ( 3 - 23 ) 

(fc,o'^(o',o) 



fc,i>0, 



Xl + A:) 2 r(l + /) 2 



We then see that in the limit (3.21) one has 

Re(Tj)^oo, i = 1,2, (3.24) 
and we can write, up to exponentially suppressed corrections, 

ziPse~ Tl , Z2~e _r2 . (3.25) 
The 't Hooft parameters of the ABJM model are large in this limit, and they behave like [7, 8] 

T1T9 

A,~^f, i = 1,2. (3.26) 

We now ask the following question: what is the behavior of the planar resolvent, or equiva- 
lently the spectral curve (3.19), in this limit? Notice that the coefficients of (3.19), regarded as an 
equation for an algebraic curve, become exponentially large or small (or they remain constant). 
This kind of behavior has been very much studied recently in the mathematical literature and 
it is known as the tropical limit of the algebraic curve (or the ultradiscretization of the algebraic 
curve), see for example [27, 28]. This limit is only non-trivial if we scale z,y in the same way, 
where 

z = log Z, y = log Y, (3.27) 
i.e. we have to consider the limit in which 



Re(z), Re(y) -»• 00 (3.28) 

as well. For generic values of z, y in this regime there is only one dominating term in (3.19), and 
the equation cannot be satisfied. To have a nontrivial equation we need at least two dominating 
terms which cancel each other. This gives us a set of linear equations on Re(z) and Re(y). 
Therefore, the "ultradiscrete" limit of the curve can be represented as a collection of segments 
in the real plane. On each of these segments there is a linear relation between z and y. It is 
an easy exercise to show that, for our particular example (3.19), the resulting diagram can be 
represented as in Fig. 2 1 . This diagram is called a tropical curve. 

This two-dimensional plane can be understood as the base of the fibration 

C* x C* -»• R 2 

(3 29) 

(Z,Y)^(\og\Z\,\og\Y\). 

1 In writing the linear equations for the segments, we have neglected constant imaginary parts, which are small 
in the tropical limit. 
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Figure 2: The strong coupling limit of the curve (3.19) can be represented as a set of segments where the 
relation between z and y is linear. This limit is nothing but the "ultradiscretization" or "tropicalization" 
of the spectral curve. 




-Ti/2«-A T x j2K.A log|Z| 



Figure 3: The tropical limit of the cuts of the Z-plane, represented here as a cylinder. 



The fiber is S 1 x S 1 and it is parametrized by the imaginary parts of (z, y). A linear relation of 
the form mz = ny + c, m,n £ Z gives a line in R 2 with a fiber S 1 C S 1 x S 1 with winding number 
(n,m). Thus the lines in the picture correspond, in the original curve, to thin tubes connected 
at the vertices. This type of picture is familiar from local mirror symmetry: as emphasized in 
[29], the mirror curve of a toric manifold, like (3.19), can be regarded as the thickening of the 
toric diagram in which lines become cylinders or tubes. In the strong coupling or large radius 
limit, the tubes become thinner and we get back the toric skeleton, which is now interpreted as 
a tropical curve. 

In order to do calculations at strong coupling we have to understand what happens to the 
period integrals in the regime (3.24). It can be shown rigorously (see for example [28]) that, in 
the limit (3.24), the periods of differentials on the original curve can be computed directly on the 
tropical curve, and they reduce to simple contour integrals along the two-dimensional diagram 
in Fig. 2. We then have to determine what is the tropical limit of the contours. To do this, we 
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Figure 4: The contour C around [1/a, a] becomes a parallelogram around the tropical curve. 



first note that in the limit (3.24) the endpoints of the cuts behave like 

A ps logC « Ti/2 ps T 2 /2, 
a ps — 6 ps £ e A , 



(3.30) 



1 1 

— 

a 6 



-A 



Let us now represent the C* domain of the variable Z as an infinite cylinder, as in Fig. 3. This 
picture also shows the contour C = C\ around the cut [1/a, a] in the tropical limit. Since our curve 
(3.19) is a double covering of this cylinder, we can build it out of two copies of C* glued along 
the cuts shown in Fig. 3. To see this in detail, let us look at the diagram shown in Fig. 2 and let 
us thicken it in order to reconstruct the original curve (3.19). If we remove the two horizontal 
segments, the thickening gives two infinite tubes which can be parametrized by z. Each of these 
tubes can be identified in turn with a copy of C*. In order to recover the full curve, we have 
to add the thickened horizontal segments. They give two horizontal tubes connecting the two 
copies of C* at z ps —A and z ps A. Notice that these locations are the positions of the small cuts 
drawn in Fig. 3. Since the solid and dashed pieces of C depicted in Fig. 3 lie on different copies 
of C*, we conclude that the contour C around the cut [1/a, a] becomes, in the tropical limit, the 
two-dimensional contour around the parallelogram shown in Fig. 4. 

Let us now consider the contour V in Fig. 1, which encircles the cut [—1/6, 1/a]. This cut 
corresponds to the horizontal tube at z ps —A, therefore the contour becomes a non-trivial cycle 
around the tube. In the tropical limit it can be schematically drawn as in Fig. 5. 

We can now use this formalism to compute some interesting physical quantities at strong 
coupling. For simplicity we will restrict ourselves to the ABJM slice N\ = N2. The resolvent is 



We first determine the relation between the 't Hooft parameter and the modulus A as follows: 



uo(z) = y(z)dz. 



(3.31) 




(3.32) 
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Figure 5: The contour V around [—1/6, 1/a] becomes a non-trivial cycle around the shrinking tube. 



which leads to 

A 71V2A (3.33) 

in agreement with the result (3.26) from [7, 8]. The 1/6 BPS Wilson loop is given by a integral 
around the contour C: 

K 8 >o=!^>- < 3 - 34 > 

In the tropical limit, this becomes an elementary integral around the parallelogram shown in 
Fig. 4. We then find, 

<"X * Si (£ - a) dz - £ < 2 + a) dz ) - £*■ (3 - 35 » 

Using (3.33) we find 

( W n 6 ) ~\V2\^ (3.36) 

which is the result obtained in [7] from the exact interpolating function, up to an overall phase 
(this is due to the fact that we neglected constant subleading imaginary pieces in the equations 
for the segments of the tropical curve). 

Another quantity that we can compute is the free energy at strong coupling. From (3.13) 
we find, 

dF 



K -7ri d> ydz = iri(p z dy « — ttiA monu y, (3.37) 
oX J v J v 

where monx> denotes the monodromy along the cycle T>. Since a ~ 6 — > 00 and Z ~ 1/a ~ 1/6, 
we have 

y « 2 log { ^JZ - 1/a - ^JZ + 1/6 J + const., (3.38) 

and 

monp y = 27ri. (3.39) 
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We conclude that 

dF 

— ^ « 2ir 2 A « 2tt 3 \/2A (3.40) 

which is the result obtained in [8]. Of course, the interest of this tropical formalism is the 
generalization to more complicated situations. This we will do in the next section, where we will 
consider the ABJM theory with fundamental matter introduced in section 2. 



4. Quenched flavor in Chern Simons— matter theories 



4.1 The quenched approximation in the matrix model 

In studying theories with fundamental matter multiplets in the context of the AdS/CFT corre- 
spondence, there have been essentially two approaches. In the first one, called the quenched or the 
probe approximation, one assumes that the number of flavors is much smaller than the number of 
colors. Since the flavor multiplets are usually obtained by adding branes to the original theory, 
the quenched approximation is equivalent to treating these branes as probes, and one assumes 
that they do not backreact on the background (see [30] for a review and a list of references for 
this approach). One can go beyond the quenched approximation and consider unquenched flavor, 
where the full backreaction of the branes is taken into account, see [31] for a recent review with 
references. It turns out that these two approaches have counterparts in the study of the matrix 
model (2.10) including flavor multiplets. We will first set up the matrix model analogue of the 
quenched approximation, and we will consider the full unquenched theory in the next section. 

In the matrix model (2.10), the inclusion of fundamental flavors leads to the insertion of two 
determinant-like operators 



(i) N 2 



n(2co Sh fp n ' 

Nf £ log (2 cosh f) - Nf £ log (2 cosh |) 



exp 



i=i 



(4.1) 



We can treat these insertions as operators which perturb the partition function without changing 
the spectral curve or resolvent of the ABJM theory. To see how this works in practice, we write 
the partition function (2.10) as a normalized vev in the ABJM theory, 



ZM=3(N 1 ,N 2 ,Nf,Nf,g s ) = (e- vv ) ABm Z A Bm(Ni,N 2 ,g s ) 



-W\ 



where 
and 



W = Nf ] Wi + Nf ] W 2 



m = ^io g 

i=i 



2 cosh ^ 



N 2 

W 2 = J^log 

3=1 



2 cosh — 
2 



(4.2) 
(4.3) 

(4.4) 



We can then calculate the free energy of the J\f = 3 theory as a cumulant expansion, 
F v - =3 (iVi,iV2,iv} 1) ,iv} 2) ,ff s ) = F AB jm(ATi, N 2 ,g s ) + ^ i_^(W fe )^ JM 



k=i 



(4.5) 
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where (c) denotes as usual the connected vev. Since W is a polynomial of degree k in the number 

(i) . . (i) 

of flavours Nj , the above cumulant expansion is an expansion around Nj = 0. Equivalently, 

we can introduce the Veneziano parameters [18] 

tf=9sNf. (4.6) 

The perturbative series (4.5) is an expansion in the Veneziano parameters around = 0, which 
is valid for 

tf < 1, (4.7) 

or equivalently 

Nf <min(iVi,iV2), (4.8) 

which corresponds indeed to a quenched approximation. Each term in this series is given by an 
integrated correlator in the ABJM theory, which is computed with the master field described by 
the resolvent (3.5). Since the spectral curve is not changed, this is equivalent to neglecting the 
backreaction of the D-branes on the original geometry. Diagramatically, the genus g correction 
to (W^abjm gives the contribution of k "quark" loops to the genus g free energy, where all 
gluon diagrams of genus g have been resummed. 

A similar perturbative scheme can be constructed for the calculation of operators O in the 
matrix model (like for example Wilson loops): 

\ / ArSJlvl ^ ^ 

The operator vevs appearing in (4.5) and (4.9) can be computed by using the connected corre- 
lation functions of the ABJM model. These are defined by 

^..->^) = (Tr^...T t ^-L_y C) (4.10) 

where U is given in (2.7). These correlators have a genus expansion 

oo 

w(z u ... ,z h ) = Y,gl g - 2+h w g {z u ... ,z h ) (4.ii) 

9=0 

which can be computed systematically with the techniques started in [32] and culminated in 
[33, 34] (for the Chern-Simons matrix models analyzed in this paper, one has to consider the 
slightly modified version of these techniques considered in [35, 36]). Let us consider for example 
the operators 

O a (X) = Trf a (X), O b (Y) = Tig b (Y), a = l,---,h 1 , b = l,---,h 2 , (4-12) 

where 

X = diag (e w ) , Y = diag (e Vi ) . (4.13) 
In this notation, the operators (4.4) are written as 

WipO =Trlog(jd , W 2 {Y) = Trlog (y* +Y~^ . (4.14) 
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We have then the following result for the connected correlators of these operators, 

(CM*) • • • o hl (x)d 1 (Y) ■ ■ ■ 6 h2 (Y))^ = 

■f 1 (X 1 )---f hl (X hl )g 1 (Y 1 )---g h2 (Y h2 ). 

This leads to a systematic 1/iV expansion by using (4.11). The planar limit of the one-point 
functions is given by the equivalent expressions 



pa p — l/b 

(O(X)) = t 1 p 1 (/x)/( Ai )d M , (O(Y)) = t 2 p 2 {v)g{v)dv. (4.16) 

Jl/a J-b 



—l/b 

pi(M)/(^)dM, (O(y)) = t 2 / 

L/ 

4.2 Quenched expansion at weak coupling 



We will now present some concrete examples of the quenched approximation, calculated at weak 
coupling. The results can be tested with perturbative calculations in the matrix model. For 

(2) 

simplicity, we will set Nj — 0, and we will focus on the free energy. The first order correction 
in Nj 1 ^ to the planar free energy is given by 

dZ 

where 



(Wi)o = - j> — log (z\ + W (Z), (4.17) 



W (Z) = ^o(Z). (4.18) 



The second order planar correction is 

and Wq(Xi,X 2 ) is the two-cut, two-point planar correlator of the matrix model. It is related to 
the Bergmann kernel of the spectral curve B(Xi,X 2 ) by [33] 

W (X 1 ,X 2 ) = B(X 1 ,X 2 ) - 1 (4.20) 

(Ai - A 2 J 

and it was first calculated by Akemann [37] in the useful form: 



Wq(XuX2) _ 1 ( l(X 1 -a)(X 1 -l/a)(X 2 + b)(X 2 + l/b) 



4(Xi - X 2 ) 2 VV (*i + b)(X 1 + 1/6)(X 2 - a)(X 2 - 1/a) 



(X 1 + 6)(Xx + l/b)(X 2 - a)(X 2 - 1/a) ^ (421) 
(Xx - a) {X 1 - 1/a) (X 2 + b) (X 2 + b) 



+ 



+ 

(a + l/b)(b + l/a) E(k) I 



4 v / < 7(X 1 )<t(X 2 ) #(fc) 2{X l -X 2 y 
where cr(Z) is given in (3.7). 
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An efficient way to calculate the above integrals at weak coupling is to perform the change 
of variables 

a — a -1 a + a -1 

X = ^— y + ^r~ ' (422) 

and expand the integrand in series in tj around ij = 0. The coefficients of the resulting series are 
relatively simple integrals, which can be computed by deforming the contour in terms of residues 
at infinity. The result one obtains is 

(Wi)o = | + ^tl (h + 6t 2 ) + ^t\t\ + g^tftj (*? - 12*1*2 + *!)+■■■, 



This can be explicitly checked against a direct calculation of the matrix integral. Indeed, we find 
in matrix model perturbation theory 



*V=3 = ^ABJM " iVj 1} 



5s 8 +5s V96" 16 192,/ 64 + 192,' + 



+ (iV 



( 1))2 



2 N[_ 3 (NfN 2 N, 
h 64 M 64 128 



, , Nf _ NfN 2 _ N?N$ _ llNf _ 3mN 2 \ 
~ \6144 ~ ~256 128 248 512 ) + 



whose planar part agrees with (4.23). 



(4.24) 



4.3 Quenched expansion at strong coupling 

Since the correlation functions (4.10) are given by contour integrals of meromorphic differentials, 
we can compute them with the tropical techniques that we introduced in the last section. We 
will focus on the free energy on S 3 , and we will assume for simplicity that iV"j 2 ^ = so that 
Tvj 1 ) = N f . We can write 

-(Wi)o = jm{z)f{z)dz (4.25) 

s c 

where 

f(z) = log (2 cosh I) . (4.26) 
In the tropical limit in which z is large the function / simplifies as 

m*Y- ( 4 - 27 ) 

Then, the contour integral (4.25) becomes 



- j> u (z)f{z)dz » J ^ {z + A) dz - J ^ (z - A) dz = A 3 « (2vr 2 A) 



3/2 (4.28) 
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and the first correction of order O(Nf) to the free energy is 



~N f Ny/2\. (4.29) 

The next order O(Nj) is much harder to compute with this technique, but it still can be 
done. Since this correction involves the two-point correlation function, which is essentially equal 
to the Bergmann kernel of the curve, what we have to do is to find the tropical limit of this 
kernel. To do this, we first discuss the tropical limit of holomorphic forms. 

In tropical geometry a tropical holomorphic 1-form is a locally constant real 1-form with 
a "conservation" condition in the vertices, and which is zero on the external legs (for basic 
notions of tropical geometry see e.g. [27]). The dimension of the space of holomorphic 1-forms 
is obviously equal to the number of independent cycles of the graph, which coincides with the 
genus of the complex curve. In our case this space is a one-dimensional space with a basis h such 
that 

h = ±dz (4.30) 
on the left and right sides of the parallelogram in Fig. 2, respectively, and 

h = ±dy (4.31) 

on the upper and lower sides. We also have h = on the external legs. 

One can realize a tropical holomorphic 1-form as a limit of a complex holomorphic 1-form: 
as we have seen, for each edge of the tropical curve we have an integer direction vector (n, m). 
Then one can associate the integral of a complex holomorphic 1-form around the corresponding 
tube with the value of the tropical one form on this vector. In this way, the "conservation" 
condition in the vertices is a trivial consequence of holomorphicity. Since the external legs of the 
graph correspond to marked points on the complex curve, the absence of poles of the complex 
holomorphic 1-form at these points corresponds to the condition that the tropical holomorphic 
1-form is zero on the external legs. In our particular case one can show explicitly that, in the 
tropical limit, 

The general notions of tropical Jacobian, Abel-Jacobi map and theta function were introduced 
in [38]. In our case the tropical Jacobian of our tropical curve C tr0 p is just 

J trop = R/LZ S 1 (4.33) 

where 



L = j> h = 8A (4.34) 
is the perimeter 2 of the parallelogram. The tropical version of the Abel-Jacobi map is 

1*trop : Ctrop ^ <Arop 

fP (4.35) 
p i—)- Utrop (p) = h mod LZ 
J Pa 

2 As usual in tropical geometry, the length of the edge is a "geometric" length with an extra weight (n 2 +m 2 ) -1 ^ 2 . 
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which equals the length of a path between the points p and pq. It can be obtained as a tropical 
limit of the ordinary Abel-Jacobi map: 



u (p) ~ ^-^trop(p)- (4.36) 



2tt 

The tropical theta function (with an odd characteristic) is 

9 trop (v) = max \ -]-Ln 2 + n(v - L / 2)\ . (4.37) 

One can easily show that 

Cp(f) = [j\ ( 4 - 38 ) 

where [•] denotes the floor function. Thus 

9" rop (v) = 6{v modLZ) . (4.39) 

This tropical theta function can be obtained as a limit of the ordinary theta function with an 
odd characteristic 

6(z) = e ^i™ 2 +27ri( 2 +r/2) _ ( 44Q ) 

In the tropical limit we have 

\L 

r« > 100, (4-41) 

2tt 



so one exponential will dominate the others in the sum. Thus one can deduce that 

' v 



lo g (i) (4.42) 



Any Bergmann kernel can be written as 

5 = B sing + 5 hol , (4.43) 

where -B s ing is given by (see e.g. [39]) 

Suing (pi,P2) = did 2 log6(u(pi) - «(p2)) (4.44) 
and i?hoi is a holomorphic part. In our case it should be chosen such that 

B = 0. (4.45) 

c 

In the tropical limit 

Bhoi(pi,P2) ~ const • h(pi)h(p 2 ) , (4.46) 

while 

£ B ing(pi,P2) ~ d l d 26'trop('"trop(P2) ~ «trop(Pl)) = -<5(«trop fe) ~ «trop(Pl))fr(Pl)fr(P2) 
= -/ldiag(Pl,P2) , 



(4.47) 
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where /idi ag is supported on the diagonal and has the property 

j> h diag (p 1 ,p 2 )f(pi,P2) = j> h(p)f(p,p). (4.48) 

CxC C 

Imposing the condition (4.45) we get 

B( P1 , P2 ) « B tIop ( Pl , P2 ) = -h diag ( Pl ,p 2 ) + h{pi) ^ P2) . (4.49) 

We can now compute the second order correction at order 0(N 2 ), (4.19), by using tropical 
techniques. It is given by 



N 2 

~2! (2vri 



' j> B( Pl ,p 2 ) log (2 cosh I) log (2 cosh I) , (4.50) 



CxC 



since the double pole subtracted in (4.20) does not contribute to the double contour integral. In 
the tropical limit this integral reads 

cic c 1 > (4.51) 

— * ^. 
3- 2 3 

Using (3.33) we obtain that the correction of order 0{N 2 ) is 

f —— A 3 / 2 . (4.52) 

96 V ; 

On the other hand, the AdS prediction for the free energy is given by (2.19). The quenched 
approximation is obtained by expanding this quantity for small Nf. Since 

^^-E (i+ Tr 3> " ^' fc -'^-f +■■■ <^> 

fe=i 

we find 

^(S 3 ) = -^N^k 1 ' 2 - ^NfNVv.+ ^NH 3 / 2 + 0(N 3 f ). (4.54) 
3 4 96 ■' 1 

We then see that the tropical computations (4.29), (4.52) reproduce correctly the first two terms 
in this expansion. 

One can try to compute the next corrections by calculating the tropical limit of the connected 
correlators (4.10) for h > 3, but as we will see in the next section it is possible to solve the planar 
theory at all values of (i.e. in the Veneziano limit) and calculate the tropical limit directly. 

5. Unquenched flavor in Chern Simons— matter theories 

We now solve the matrix model (2.10) in the planar limit, but for all values of by using the 
techniques of [10, 14]. 
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5.1 Exact resolvent in the Veneziano limit 

The starting point in the calculation of the resolvent are the saddle-point equations 
* + \- tanh m=T, coth - V tanh , 

r(2) 7V 2 JVi 



(5.1) 



— - + 



f — tanh (^) = j2 coth ^ 2 ~ " S tanh 



q s 2 V 2 / ^ 2 ^ 2 

ys fe^a i=l 

We will solve instead the problem 



^ + !i_ tanh g)4r C oth + tanh 

p 2 V 2 / JVi ^ 2 N 2 ^ 2 

a=l 



Hi - 



f a — tanh ( — ] =-J- coth — — + —J- tanh — — . 

2 V 2 / iV 2 ^ 2 iVi ^ 2 



(5.2) 



(1 2) 

analytically in the parameters ti j2 and tf , and then we will perform the analytic continuation 

t 2 -»• -t 2 . (5.3) 

The procedure to solve this type of equations is as in [9, 10, 14]. We first introduce exponentiated 
variables 

Zi = ew W a = er°. (5.4) 
In terms of these variables the saddle-point equations read 

W ry nj . iVl 0i AT 2 



ty_Zi-l JVi-l 2ti^ Z, 2t 2 ^ Wq 

J^l J o=l 

if Wo - 1 , iV 2 - 1 2t 2 ^ w b 2*i A Zi 



(5.5) 



. V VK a - 1 , , , iV 2 - 1 Zl 2 Wf> ^1 



2 W a + 1 -/V 2 iV 2 ^ W a - W b Nx^Wa + Zi 

byta i=l 



The resolvent ooq{Z) is defined as in (3.4), and it will have two cuts corresponding to the set of 
eigenvalues. Let [a, b], [c, d] C M be the cuts corresponding to — W a and Zj, respectively. In terms 
of the planar resolvent we have, 



t w z _ 1 1 

logZ+^-^ I =-( Wo (Z + iO) + Wo (Z-iO)), 
log(-W) - ^-^^ =2 + iO) + wo(W - i0)) 



(5.6) 



As in [14] we now define the functions, 

F(Z) = yfi(ZjJ*dX 



G(Z) = y/<r(Z) / dX 

J a 



Z - X' 

g(x) 

Z -X' 



(5.7) 
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where 



fix) = . 1 \\ogx + \ 

(5.8) 

^ 1 ^^ ( ^ t f x + l 
oi = . log(— x) 

It is then easy to show that the planar resolvent, defined again by the VEV (3.4), is given by 

u (Z) = -(F(Z) + G(Z)). (5.9) 
One can check that, as in the ABJM matrix model, 

ab = l, cd=l. (5.10) 

This follows from the symmetry of the saddle point equations under Z — > Z -1 , W — > W^ 1 , 
together with the conditions on the endpoints imposed by the asymptotic behavior of ujq(z) at 
infinity, 

loo(z) ~ t, z — y oo. (5-11) 

From now on we will denote the two independent endpoints as a, 6, as in the ABJM model, so 
that the cuts are [1/a, a] and [—6, —1/6]. The two equations 

r dZ u (z) i dZu^z) 

h = f Cl ^^r> t2 = L^— (5 - 12) 

determine them a, b as a function of ti,<2- Here, Cj encircle [1/a, a] and [—6,-1/6], respectively. 
We can now calculate the planar resolvent explicitly. It is given by two pieces. The first one 

is, 

dX log A f' 1 /" dX log(-X) 

(5.13) 



ti/a z - X y/\^X)\ J. b z-x^xJi 1 
which is simply the resolvent of the lens space matrix model (3.5). The second piece is 

dX X-l 1 tf r-i/b dX X + 1 1 



t)' r a dx x-i i ty r 
~2~ Lz-ii + iJ^w + T J 



/a Z-XX + l^(x)\ 2 J_ b Z-XX-l^(x)\ 



(5.14) 



These integrals can be expressed in terms of elliptic functions. In order to do so we use the 
results 



/_ 



A/a (£- A) v / |cr(A)| l + a6(Z-a)(Z + 6) 

<L ' ; (-(a + 6)n(n-(Z),fe) + (Z + 6)^(fc)) 



/-& (Z - X)yJ\a{X)\ l + ab{Z - a){Z + 6) 

(5.15) 

where 

2 _ (a 2 -l)(6 2 -l) n+(zl _ a 2 ~l^ + 6 () _ _6 2 _-l_Z-a 

fc " (l + a&)2 ' ( j " l + a6Z-a' " [Z) " l + a6Z + 6' ^ 
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Defining the auxiliary function 
J(a, b, Z, s) = 

2y, " h ' ^(a + b) - Z ~° ■ mn+(Z),k) + ^K(k))-(Z^-.^ ^ 



18) 



l + abZ + s\\ y J (Z-a){Z + b) y v n ' Z + b 
we finally obtain 



uo(Z) = cj£ bjm (Z) + ^5 J(a> 6) Zj i) + *(?) J( _ 6j _ 0> z> ; (5 . 



where Wo lBJM (Z) is the resolvent in the theory without matter, and it is given in (3.5). The 
asymptotic behavior (5.11) determines 

t = log (/?) 

- ^Y^ab) ( f ? ((6 + l)K{k) ~ {a + 6)II(na ' k)) + (_(a + l)K{k) + (a + 6)IT(n6 ' ' 

(5.19) 

where we have used the notation 

I -a 2 1-b 2 , , 

na = rr^' nfe = rr^- (5 - 20) 

The relation (5.19) reduces to (3.12) when both t^ and t^ go to zero. 
5.2 Weak coupling limit in the unquenched theory 

In order to test the above expressions, we can compute the expansion of the resolvent at weak 
't Hooft coupling (i.e. around ti = 0) but for arbitrary To do this, the first step is to express 
the endpoints of the cuts in terms of the 't Hooft parameters. The period integrals for ti can 
be expanded around a = 1, b = 1, and these series expansions can be inverted. At the first few 
orders in ti we find 




1 if 2fa+ 6 (2f) 8/a v t h ^ /2+ (Tf) 3/2 ^ 




(5.21) 



+ , \— 7 I — ! + — — - i-i + 
2 (if ) 

In these equations, 



^ = 1 + ^-, Tj 2) = l-^-. (5.22) 
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When tf = we recover the mirror map at the orbifold point of [9]. As a test of these results, 
we can calculate the coefficient of 1/Z in the resolvent. After the analytic continuation (5.3), 
this coefficient computes the planar VEV of the supertrace of U in the matrix model (2.10), as 
in (2.6): 

2g s (Str D U) = ( 
t {1) 

+ ^ 7 ((! + ab) 2 E{k) + {b-a + a 2 b + 3ab 2 )K(k) - 4a6(a + b)U(n a , k)) 

27rV ab(\ + ab) (5.23) 

t m 

+ jJ- ((1 + ab) 2 E{k) + (a-b + ab 2 + 3a 2 b)K(k) - 4ab(a + b)U(n b , k)) , 

27rva&(l + ab) 

where (5.3) must be implemented at the end of the calculation. At weak coupling we find, 



g s {Sir n U) Q = h + t 2 H 



2Ty> 2T f 

1 I 1 | 3 t\t 2 htl 1 / 1 ! 

(5.24) 

The r.h.s. has the expected symmetries of the matrix model. Indeed, it is odd under 

h o -t 2 , tf o -tf. (5.25) 

Notice that each term in the expansion (5.24) is a rational function of the Veneziano parameters. 
From the diagrammatic point of view, each of these terms corresponds to a fixed planar "gluon" 
diagram (with a boundary associated to the insertion of U) in which we have summed over all the 
"quark" loops, i.e. the "gluons" are quenched and the "quarks" are dynamical. One can actually 
check the first few terms written down in (5.24) againts an explicit perturbative calculation in 
the matrix model. 

Based on [4], one should expect that the VEV (5.24) computes (twice) the VEV of the 1/2 
BPS Wilson loop operator. However, in order to assert this one should first check that the 
construction of [4] of this operator extends to the flavored theory that we are considering here. 
In any case, the formulae (3.13), (3.14) remain valid in the flavored theory, since the 1/6 BPS 
Wilson loop operator can be constructed for all Chern-Simons-matter theories considered in 
[2]. We will now evaluate these formulae in the unquenched theory at strong coupling, by using 
tropical techniques. 

5.3 Strong coupling limit in the unquenched theory 

For simplicity we will set N^f = 0, Nj = Nf- We will write (5.18) as 

uo(z) = y(z)dz, y{z) = y p (z) + fiy m (z). (5.26) 
In this equation, /j, is defined in (2.17), 

y p (z)= Wo ABJM (z) (5.27) 
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A-K 



-(A-K) 




-A o A 




Figure 6: The two-dimensional graphs representing the tropical limits of y p (z) (left) and y m (z) (right). 



is the equation of the spectral curve (3.5) in the ABJM model, and 

, , 1 fdX 1 X - 1 JaiZ) 
Vm{z) ~ 1 



2 J 2vri Z - X X + 1 y/^(Xj ' 



In the tropical limit £ pa e A is large, and we will set 



K 



(5.28) 



(5.29) 



where K is a parameter to be determined. In the ABJM model with N\ = N2 one has K = 0. 
We will shortly determine the value of K in the theory with unquenched flavor, in the tropical 
limit. We will assume that A > 0, \K\ < A, which will be justified a posteriori. 

The integral (5.28) can be evaluated "tropically." To do this we use the tropical limit (4.32) 
of the holomorphic one-form, as well as the limits 



Z J 1, ifx<z 

Z - X - e {x z > , if x > z 



and 



±E(z) = ± { 



1, if -A < z < ^4 



(5.30) 



(5.31) 



^z-A 



-e" if z > v4 

Here the ± sign corresponds to the two determinations of the square root. One finally obtains 



y m (z) « ± ^ H(z) ^ /i x $(x, z) signx, 



(5.32) 



where fo x is defined in (4.32) but with z replaced by x. In this result, signx/2 can be interpreted 
as the derivative of \x\/2, which is the tropical limit of the potential deformation, see (4.27), 
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Figure 7: The tropical curve representing the tropical limit of the resolvent (5.18) of the Af = 3 theory. 

and in principle we can generalized it to other deformations. The tropical limit of y m (z) can be 
rewritten as 

, s , 1 /, . (A-K+(A-\z\), if-A<z<A, 

y m {z) ~ ± - n x signx = ± < (5. do) 

2 J \ A — K, otherwise, 

c z 

where C z is a line connecting two different points on the curve with the same value of z. Equiv- 
alently one can obtain (5.33) by taking the tropical limit in the explicit expression (5.18). The 
tropical limit of the two-valued functions y P (z), ym(z) can be represented by the two-dimensional 
graphs shown in Fig. 6. Of course, the diagram for y v {z) is nothing but the tropical curve rep- 
resented in Fig. 2. 

We now want to find the relation between A and K. To do this, we will impose for simplicity 
that N\ = N2 = N in the M = 3 Chern-Simons-matter theory. This means that the total 't 
Hooft parameter t = t\ + ti vanishes. It follows from (5.12) that this sum can be evaluated by 
deforming the sum of the contours C\ and Ci to infinity and the origin, so we obtain 

resz=oo wq(Z) - res z =o uq(Z) = 0, (5.34) 

which leads to 

K = -^A. (5.35) 

Notice that in the limit /j = 0we correctly reproduce K = 0. With this relation at hand we can 
already add the two graphs to obtain the tropical curve representing the new resolvent ujq(z), 
which is shown in Fig. 7. The calculation of the different periods reduces, like before, to trivial 
line integrals on the plane. We first have to relate A to the 't Hooft parameter. We have, for the 
period (5.12), 

1 / , s (1 + M/2M 2 
2vriA = ti = —<b y{z) dz ps --^ ^—^ — (5.36) 

47T1 J c 7T1 
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and we find 

A ^ , (5.37) 

which is the deformation of the relationship (3.33) in the ABJM theory. For the planar free 
energy we have, as in the period integral (3.37) for the ABJM theory, 

= -iriAmonvy. (5.38) 

The monodromy of y m (z) can be computed as (3.38). Indeed, the main contribution to the 
integral for y m is given by the vicinity of 1/a, and we have 

dX y/(Z-l/a)(Z + l/b) 



< a X-Zy/V-V*){X + W (5 .3 9 ) 
= log | ^Z - 1/a - y/Z + 1/6} - log [y/Z - 1/a + yj Z + . 

Then the monodromy around the cut [—1/6, 1/a] is 

monx)y m ~ 27ri . (5.40) 

We conclude that 



or equivalently, 



2ttM(1 + /x), (5.41) 



F (A)^vr 3 ^ 1±^ . (5.42) 

This is in perfect agreement with the AdS prediction (2.19). Notice in particular that we have 
been able to reconstruct the full nontrivial function £(//) involved in the volume of the tri-Sasakian 
target space (2.16). 

Finally, we can calculate the vev of the 1/6 BPS Wilson loop, which is given again by (3.34) 
but now with the new resolvent. We obtain, 



W r 



° /6 )o ^(f; 2 ^-^^- N)) dz-j\*(z + A + - \z\)) dz) 

(5.43) 



2A / / 2A 

exp 7T4 



2 y i + n/2 \ y i + n/2 1 

With some more work one can show that the would-be 1/2 BPS Wilson loop has the same leading, 
exponential dependence. This is in perfect agreement with the AdS prediction (2.21), and the 
exponent should be equal to the regularized area of a fundamental string in the corresponding 
type IIA background. 
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6. Conclusions 

The two main results of this paper are the following. First, we have shown that the strong 
coupling limit of Chern-Simons-matter theories corresponds, in the matrix models of [2], to the 
tropical limit of the spectral curve. In this limit, the curve becomes a two-dimensional graph 
and the period integrals which give the main observables of the matrix model can be computed 
in a very simple way by evaluating line integrals along the graph. Second, we have obtained 
the planar resolvent of the matrix model describing the ABJM theories with flavor introduced 
in [15, 16, 17], including all planar loops of fundamental matter. In other words, we have solved 
the Veneziano limit of these flavored theories. The unquenched regime at strong coupling can be 
studied by applying our tropical techniques to the exact resolvent, and we have reproduced in 
detail the AdS predictions of the conjectural large N duals. In particular, we have reconstructed 
from a gauge theory calculation the volume of the tri-Sasaki Einstein manifolds involved in these 
AdS duals. 

There are many interesting problems which deserve further research. It would be interesting 
to analyze with our tropical techniques other Chern-Simons-matter theories, like for example the 
two-node quiver with arbitrary levels introduced and studied in [40]. The exact planar resolvent 
was recently written down in [14], and our techniques should lead to explicit expressions for its 
free energy at strong coupling. Since this theory has a conjectural large iV dual described by 
massive type IIA supergravity, one might be able to test directly in the matrix model many of 
the results recently discussed in [41]. 

On a more mathematical level, we have seen that the quenched approximation to the flavored 
theories studied in this paper involves in a natural way the tropical limit of the connected matrix 
model correlators (4.10). These correlators can be computed with the topological recursion of 
Eynard and Orantin [34], and in section 4 we have described the tropical limit of one of the 
basic building blocks of this recursion, namely the tropical Bergmann kernel. It would be very 
interesting to see if the full topological recursion can be formulated directly in the tropical limit. 
This would give a powerful formalism to study perturbations of ABJM theory and it might find 
other applications in the future. 

It seems that, in the Chern-Simons-matter theories with less supersymmetry, the planar 
resolvent is much more complicated than in ABJM theory, and in particular we don't longer have 
the underlying local Calabi-Yau geometry that we had in (3.19). For example, the resolvent in 
(5.18) involves elliptic functions, in contrast to (3.5). The tropical limit is however given by a 
graph which ressembles very much a blown- up toric Calabi-Yau manifold (see Fig. 7). It would 
be important to understand at which extent the solution to these Cherm-Simons-matter matrix 
models can be encoded in a toric Calabi-Yau. 
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A. Strongly coupled density of eigenvalues and tropical geometry 



In this Appendix we rederive some of the results for the N = 3 theory by using the approach of 
[19], and we compare it in detail to our tropical methods. 

The starting point of [19] is an analysis of the ABJM matrix model (2.1) in the ABJM slice, 
at large N but fixed k, which corresponds to the strongly coupled limit of the theory. This 
analysis makes possible to re-derive in a particularly simple and beautiful way the result of [8]. 
Let us see how this is done, following closely the steps in [19] . The behaviour at large N of the 
equilibrium eigenvalues of the matrix model is, as derived in [8], 

p k = N 1 / 2 x k + i£ k , v k = N 1 l 2 x k -i£ k , k = l,---,N, (A.l) 

where x k , l k are of order one at large N. At large N the eigenvalues x k , i k become dense, so 
that 

^ -> i e [0, 1] (A.2) 

and they are described by the functions 

p(x) = ^, £(x). (A.3) 
It is shown in [19] that, when N is large, the free energy of the matrix model can be written as 



iV 3 / 2 



— f dx xp(x)£(x) + f dx p 2 {x)f {2£{x)) — — [ / dx p(.r) 



(A.4) 



Here, f(t) is a periodic function of t, with period 2ir, and given by 

f(t) = ir 2 -t 2 , i€[-7r,7r]. (A.5) 

The last term in (A.4) involves, as usual, a Lagrange multiplier m imposing the normalization of 
p(x). Notice that our sign convention for the free energy is the same as in [8] and opposite to the 
one chosen in [19]. Varying this functional w.r.t. p(x) and £(x) one obtains the two equations 

2irp(x)f'(2£(x)) = -kx, 
4vrp(x)/ (2£{x)) = m - 2kx£(x), 

which are solved by 

= <w = ^ (A - 7) 

The support of p(x), £(x) is the interval [— x*,x*]. One fixes and m from the normalization 
of p and by minimizing — F. This gives 

/2 2vr 3 , 

x * =7r \k 7 m= ^7' ( ^ 

Evaluating the free energy for the functions (A. 7) and the values (A. 8) of x*, m, one reproduces 
the result of [8] for the free energy. 

The above results can be easily compared with our tropical analysis. The value of x* gives 
(up to a factor A 1 / 2 ) the position of the endpoint A, and it is in accord with the value of (3.33), 
since 

A = N l > 2 x*. (A.9) 
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The fact that the density p(x) is constant follows from our result for the tropical limit of the 
curve. Indeed, the density of eigenvalues (normalized as to have an integral along the cut equal 
to one) is given by the well-known formula 

p(z) = ^disc y(z) (A.10) 

where discy(z) is the discontinuity of the curve through the cut [—A, A]. In our case this is just 
the constant 2A, and it is given by the horizontal separation between the two diagonals in Fig. 4. 
Changing variables from z = N 1 / 2 x to x we find indeed, 



\/2k 

P ( X ) = (A.11) 

in precise agreement with the result of [19]. 

The inclusion of fundamental matter in the approach of [19] is straightforward. — F includes 
now the extra term 

J dxp(x)\x\ (A.12) 

which is the large N limit of the operator in the exponential of (4.1) (as in (4.27)). The new 
saddle point equations are 



2np(x)f (2£(x)) = -kx, 
A-np(x)f {2i{x)) = m - 2kxt{x) - *Nf\x\, 

with solution 

m — irNf\x\ „, s kir 2 x 



(A.13) 



p{ x ) = /' £( x ) = — _ _ . (A.14) 

Normalization of the density and minimization of — F lead to 

z. = —*^ , m = 2^4±£^. (A.15) 
^2k + Nj ^/2k + Nj 

A straightforward calculation of —F reproduces (2.19). 

Let us now compare this with the tropical approach. First of all, we have again the equality 
(A. 9) between the endpoints of the cut in both approaches, involving now the value of x* obtained 
in (A.15) and the value of A obtained in (5.37). The density of eigenvalues in the theory with 
matter can be obtained from the planar resolvent (5.26) as 

ATl/2 

P(x) = ^2^ dlSC (Vp( Z ) + Wm(^)) , (A.16) 

where p is defined in (2.17). From Fig. 6 and (5.33) we read off immediately 

^discy p (z) = A-K, -discy m (z) = 2A-K-\z\, (A.17) 

and we deduce 

P(x) = ^ [(1 + P)A - pN^M] ■ (A.18) 
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Plugging in the value of A (5.37), we recover precisely the form of p{x) given in (A. 14). This 
shows explicitly that the piece- wise linear densities obtained with the method in [19] correspond 
to the tropical curves obtained in this paper. 

One advantage of the method of [19] is that it gives the large N limit of the free energy, 
at strong coupling, without directly using the resolvent of the model. Therefore, this method 
is useful when the resolvent is difficult to write down. If on the contrary one is interested in 
calculating the resolvent of the model (to study for example weak coupling expansions), then 
the tropical approach developed in this paper provides a powerful method to extract the strong 
coupling limit from the resolvent. In addition, the method of [19] assumes as an ansatz the scaling 
(A.l) of the eigenvalues, while the approach based on studying the resolvent and its tropical limit 
provides a bona fide solution to the large N theory without any further assumptions. 
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